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Abstract 

Once a light Higgs boson is discovered at a future e + e~ collider, the next target at 
the collider will be precise measurements of the Higgs boson properties. In this paper 
we report a simulation study on the measurement of the ratio Br(H — ► cc, gg)/Br(H — ► 
bb) at center of mass energy of 300 GeV, and show the possibility to constrain MSSM 
parameters from the measurement. 
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1 Introduction 



Next generation e + e~ linear colliders, firstly to be operated at center of mass energies (a/s) 
of 300-500 GeV (phase 1) and to be then upgraded to 1 TeV or higher (phase 2), have been 
proposed by several groups in recent years [[TJ 0, f|, [50. At the first phase of such colliders 
the main target is the Higgs boson, the only undiscovered participant in the Standard Model 
(SM). According to the recent precision electroweak measurements at LEP and SLC, Super- 
symmetry (SUSY) has become more favorable. In the Minimal Supersymmetric extension of 
the Standard Model (MSSM), two Higgs doublets are necessary to give the masses of up and 
down type fermions. Two neutral scalars (h°, H°), a pseudoscalar (A ) and a pair of charged 
Higgs (H^) are introduced in this extended Higgs sector. The mass of the lightest neutral 
Higgs (m h o) is smaller than that of the Z° boson (m z o) at tree level. Even if radiative correc- 
tions are taken into account, m h o is at most ~ 130 GeV for a top mass (m t ) of 170 GeV 0, ||. 
This situation does not change significantly in extended SUSY models [^]. Therefore if SUSY 
is correct, discovery of at least one neutral Higgs boson is guaranteed at the first phase. 

After having discovered the Higgs bosonQ, the next step is to distinguish whether the 
Higgs is SM-like or SUSY-like. This can be done by measuring the ratio of branching frac- 
tions Br(H — > cc)/Br(H — > 66) or almost equivalently Br(H — > cc, gg) / Br(H — > 66)0. In the 
MSSM these values are smaller than those of the SM. As the value Br(H — > cc)/Br(H — > bb) 
(or Br[H — > cc, gg)/Br(H — > 66)) has a strong correlation to the A mass (m A o) and al- 
most independent from the SUSY energy scale (m SU sY)/ it is pointed out that a constraint 
on m^o can be obtained from this measurement [[10]]. Fig. p] shows the relation between 
Br(H — > cc, gg)/ Br(H — > bb) and m^o. 
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Figure 1: Relation between Br(H —>■ cc, gg)/ Br(H — * bb) and m^o [1C]. 



The mass of the bottom quark (m&) is important in this relation, since the ratio Br(H — > 
cc, gg)/Br(H — > 66) is proportional to m^ 2 . We can use the measurement of Br(H — > t + t~) 
to determine ra&, through the relation Br(H — > t + t~) / Br(H — > 66) oc m^/mf which is valid 
both in the SM and SUSY models. 

1 Hereafter we denote H as the discovered Higgs boson. 
2 We define Br(H — > cc, gg) as Br(H — > cc) + Br{H — > 53). 
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The total production cross section of the process e + e — > HZ (a tot ) can be determined 
by counting the events, in which the Z° decays into a pair of leptons whose recoil mass is 
consistent with the Higgs mass. This technique is valid even if H has decay channels to 
undetectable particles such as neutralinos, since only the Z° decay products are used in the 
event selection. 

A previous simulation study for the measurement of the Higgs boson branching frac- 
tions can be found in Ref . where the measurement was studied at y/s = 400-500 GeV 
mainly for tuh = 140 GeV. The authors of Ref. [|11|] claimed that the statistical error of the 
Br(H — > cc, gg)/ Br(H — > bb) measurement would be about 40 % for the 120 GeV Higgs 
with an integrated luminosity of 50 fb ~ 1 . However, their result is not sufficient to provide 
strong constraints on m^o, and the Higgs mass of 140 GeV may be too heavy if we assume 
the MSSM. Therefore we perform here an extended simulation study focusing on the mea- 
surement of the ratio Br(H — > cc, gg)/Br(H — > bb) assuming a SM-like Higgs, where we try 
to improve the measurement as follows; 

• setting y/s =300 GeV to avoid large background from top pair productions 

• using all the Z° decay modes to reduce statistical errors 

• studying the effect of beam polarization to reduce electroweak background. 



2 Event Generation and Detector Simulation 



We use PYTHIA5.7 + JETSET7.4 [|12|] to calculate the cross sections and to generate events of 
the processes e + e~ — > Z°H and major background processes. We use the quark masses of 
m t = 170 GeV, m b = 4.25 GeV and m c = 1.35 GeV and initial state radiations are considered. 
A cross check using GRACE 1.1 [ |T3[ ] shows no inconsistency in cross sections and angular 
distributions. Table [1] shows the list of branching ratios, cross sections and expected number 
of events in a year (50 fb -1 ) for the SM Higgs boson. Number of generated events both for 
signal and background processes in this study are also listed. The cross sections with 90 
% polarized electron beamf] are obtained from the cross sections and angular distributions 
calculated with GRACE 1.1. 

The 4-momenta of all decay products given by the event generator are fed into a fast 
detector simulation program which is based on a smearing method with the standard pa- 
rameters of the JLC-I Detector given in Ref. [|IJ]. The characteristics of the detector simulator 
are the following: 

• A pixel type vertex detector (VTX) with a pixel size of 25 x 25/oti 2 is assumed. 

• The momentum resolution of the central drift chamber (CDC) is assumed to be cr pt /pt = 

1.1 x 10~ 4 • p t © 0.1% (p t in GeV). 

• The electromagnetic and hadronic calorimeter hits are created according to the energy 
resolutions of 15%/ VE © 1% and 40%/ VE © 2%, respectively (E in GeV). 

• The momentum measured with the CDC is always used for a charged track rather than 
the energy measured with the calorimeter. 



3 Polarization of electron is defined as: P e - = (N^ — /(N^ + N-&), where JVl and TVr are the number of 
left-handed and right-handed electrons, respectively. 
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• The contribution of charged tracks is subtracted from the energy of the calorimeter 
clusters in order to extract the energy of neutral particles. 

• All electrons (muons) with momentum above 1 GeV(2 GeV) are assumed to be perfectly 
identified. 



Table 1: Branching ratios, cross sections, expected numbers of events, expected numbers of 
events with 90 % polarized electron beam and numbers of generated signal and background 
events in a year for the SM Higgs boson with a mass of 120 GeV (50 fb _1 ). Cross sections 
for e+e - — > euW, e + e~Z° with 90 % polarized electron beam are not calculated here since no 
events survived the selection criteria described in section [|. 





120 GeV 


Br (%) 


a(fb)xBr 


events /yr 


90% Pol. 


Created 


e+e" 


-> Z"H 




191 


9,540 


7,920 




H 


-> bb 


66.2 


126 


6,310 


5,240 


50,000 


H 


— > cc 


2.40 


4.58 


429 


190 


10,000 


H 


99 


8.33 


15.9 


794 


659 


10,000 


H 


-> cc + gg 


10.7 


20.5 


1020 


850 


10,000 


H 


— > T + T~ 


7.28 


13.9 


694 


576 


10,000 


H 


-> w + w~ 


13.4 


25.6 


1,280 


1060 


10,000 



e+e" -> W+W- — 13,900 697,000 80,200 700,000 

e+e" -> Z°Z Q — 1,080 53,900 36,100 108,000 

e+e" -> e+e-Z° — 8,310 416,000 — 100,000 

e+e" -> ez/H/ — 2,290 115,000 — 100,000 



3 Higgs Event Selection 

The final states of the process e+e - — > can be classified into three types in terms of the 
Z° decay modes. 

1) 4 jet Z°^qq -70% 

2) 2 jet Z° -> i/p ~ 20 % 

3) 2 jet+ZZ Z°^l+l~ -10% 

The event selection for each final state is described in this section. The selection criteria 



for the a tot determination is described in section p3| . All cut values described here are for 
the Higgs mass of 120 GeV. Some of the cut values are changed for other Higgs masses to 
optimize the signal/background ratio. 

3.1 4 Jet Analysis 

In the 4 jet analysis, each event is forced to be reconstructed as a 4 jet final state (two jets 
from Higgs and other two from Z°) with the JADE jet finder []T4l]. All possible combinations 



of jets are examined to select the combination which minimize the value 

(Mij - m H f + (M kl - m z o) 2 i,j, k,l = l, 2, 3, 4, 
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Figure 2: Distribution of Thrust value for the signal and background events in 4 
jet analysis. 



where M^- is the invariant mass reconstructed from the jets i,j assigned to the Higgs jets, 
and M ki is that from the jets k, I assigned to the Z° jets. We denote M& and Mjf t as the 
reconstructed Higgs and Z° masses in the selected combination. Cuts on the visible energy 
of E vis > 0.85 • y/s and on the longitudinal and transverse momentum balances of | J2Pz I < 30 
GeV and \J2Pt\ < 30 GeV are then applied. A containment cut of | cos thrust I < 0.8 and a cut 
on the thrust value Thrust < 0.75 are applied to reduce e + e _ — > W + W~ , Z°Z° backgrounds 
(see Fig. 0). Events in which the value A Z h defined as 



is less than 5 are selected as Z°H candidates. Once events are selected as Z°H events, we ap- 
ply further cuts to select signal events (H — > bb, cc and gg). To reduce H — > t + t~ background 
events, the number of charged tracks from the Higgs jets is required to be greater than 20. 
This rather large cut value is chosen because there is inevitable contamination of tracks from 
the Z° jets to the H jets in the jet finding procedure. The value returned from the JADE 
jet finder for 6 jet final state, must be smaller than 0.005 to reduce the H — > W + W~ back- 
ground. This cut is based on the fact that the background event H — > W + W~ has 6 jets in 
the final state while the signal events have only 4 jets. The expected numbers of events to be 
collected in a year and the selection efficiencies^ are summarized in Table 0. 




4 Selection efficiency is defined as: e sc i 



# of selected events 



# of generated events 
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3.2 2 Jet Analysis 



In the case where the Z° decays to vv, the final state has two jets from the Higgs accompanied 
by a large missing energy Since the reaction e + e~ — > Z°H is essentially a two body process, 
the missing energy distribution has a sharp peak although initial state radiations degrade 
this character. The events are required to have a missing energy of 139 GeV < E mis < 149 
GeV (see Fig.^j). A containment cut of | cos 0j et | <0.8 is also applied. Events with high energy 
leptons (-Eiepton > 10 GeV) are removed. The event is accepted as a Z°H candidate if the 
value A Z h defined as 



A 



ZH 



jet ~ m H 



+ 



^r 1 



(J n,f recoil 
M jet 



is less than 4.5, where M- ]et and Mj r e c t C011 are the invariant mass and the recoil mass recon- 
structed from the two jets. After selecting the Z°H candidates, cuts of N track > 10 and 



Ycut < 0.012 are applied to reduce the H 



, W + W backgrounds, respectively, simi- 



lar to the case of the 4 jet analysis. The expected numbers of events in a year and the selection 
efficiencies are summarized in Table 0. 
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Figure 3: Missing energy distribution for the signal and background events for 2 
jet analysis. 



3.3 2 Jet + // Analysis 

In this analysis, two same flavor leptons (e, jj) with opposite signs are required in an event. 
After removing the two lepton tracks, all events are reconstructed as a 2 jet final state. Con- 
tainment cuts of | cosset | <0.8 and | cos lepton | <0.8 are applied. Events which satisfy the 
condition 

Xzif= j Mf^n - rriH V | / M lepton - m z o \ 2 < ? 
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Figure 4: M^"^ l n distribution after all cuts for 2 jet+ZZ analysis. The shaded his- 
tograms are for the background events. 



are selected as Z°H candidates, where M lepton and M^^ n are the invariant mass and the re- 
coil mass reconstructed from the lepton pair. We use here M{^ l n as the Higgs mass. The 
same cuts as the 2 jet analysis are applied to reduce the background events from H — > 
t + t~, W + W~ . The distribution of in this analysis is shown in Fig. f|. 

The selection criteria for the a to t measurement is almost the same as that of 2jet+ZZ analysis 
except that the cut for cos 9 ]Ct is not used. The expected numbers of events in a year and the 
selection efficiencies are summarized in Table @. 

Table 2: The event selection efficiencies and expected numbers of the Z°H events in a year 
(50 fb _1 ) for the SM Higgs boson with a mass of 120 GeV. Those of the background processes 
are also given. 



4 jet 2 jet 2jet+/Z a tot 

e sc \ (%) # events e sc \ (%) # events e sc \ (%) # events e sc \ (%) # events 

H 

H 

H 

H 

H 

H_ 
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3.0 


99 


12.8 


102 


4.5 


35.4 
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0.0 
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0.0 
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70.5 


e + e'Z° 
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0.0 


evW 


0.0 
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0.0 


0.0 


0.0 


0.0 


0.0 


0.0 
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4 Flavor Tagging 



The technique of the flavor tagging used here is based on the long lifetimes of hadrons 
containing heavy flavor quarks. Tracks from the decay of such hadrons have large impact 
parameters (b) from the primary vertex. The impact parameter can be precisely measured 
with the VTX surrounding the beam pipe. 

In this study, we calculate the three dimensional impact parameter using the space point 
information from the VTX and the momentum information from the CDC. The error of the 
impact parameter measurement (a b ) is calculated according to the following formula [[!]]. 



\ 



11.4 2 + 



28.8 



^pV sin 3 9 , 



(/jm) (p in GeV) 



The first term in the square root is the geometrical contribution which is determined by 
the pixel size and the geometry of the VTX, while the second term is the contribution of 
multiple scattering which depends on the thickness of the detector. We use the normalized 
impact parameter b/a b to scale the impact parameter, instead of the impact parameter itself. 
The number of tracks in the Higgs decay, with their normalized impact parameters being 
greater than a cut , is used for the flavor tagging. Fig. [5| shows the distribution of the number 
of such tracks for the three Higgs decay modes together with the selection regions. The 
value a cut and the selection regions are optimized to minimize the statistical error of Br[H — > 
cc, gg)/ Br(H — > bb). We use a cut = 3 here. The tagging efficiency^ and the expected number 
of tagged events in a year are summarized in Table ^ for each final state. 
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Figure 5: Numbers of large impact parameter tracks in jets from the Higgs decay 
for three Higgs decay modes (A): H — > bb, (B): H —> cc and (C): H — > gg. 
The outside events of solid line are assigned to bb, and inside ones are 
assigned to cc, gg. Inside of the dashed line is assigned to cc. 



^Tagging efficiency is defined as: e ta 



# of tagged events 
# of Z°H events 
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Table 3: The tagging efficiencies and expected numbers of signal and background events for 
the SM Higgs boson with a mass of 120 GeV (50 fb _1 ). Here a cut = 3 are used. Efficiencies 
are given in percent. The bold numbers are efficiencies of the signal events. 
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5 Selection of H — ► r + r 

The event of the process e + e — > with the decay H — ► r + r has at least two unde- 
tectable neutrinos in the final state. Therefore the Z°H event of this decay mode should be 
identified using mainly Z° decay product, and the 2 jet events where Z° decays into two 
neutrinos are not used in this analysis. The strategy of the event selection is the following: 

• Select Z°H event using the Z° decay product and the opening angle (8 H ) of the jets 
from the Higgs decay. 

• Select events with a decay H — > t + t~ by using the missing energy and the charged 
multiplicity. 

In the 4 jet analysis, we select a pair of jets whose invariant mass is closest to m z o . Cuts 
of | cos thrust I < 0.8, Thrust < 0.8, E mis > 10 GeV, and 9 H > 90° are then applied. We require 
that the recoil mass against the Z° is consistent with the Higgs mass, and that each jet from 
the Higgs decay contains at most five charged tracks. 

In the 2 jet+ZZ analysis, we require that there is a pair of opposite sign leptons in the event, 
and that the invariant mass and the recoil mass of the leptons are consistent with Z° and H, 
respectively. Cut on the missing energy of E mis > 10 GeV is then applied. The two jets 
are considered to be the Higgs decay product and the number of charged tracks in each jet 
should be less than three. The selection efficiencies and the expected numbers of signal and 
background events in a year are summarized in Table f|. 
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Table 4: The selection efficiencies and expected numbers of the H — > t + t events for the SM 
Higgs boson with a mass of 120 GeV (50 fb _1 ). Those of background events are also given. 
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6 Results 

We list in Tables @, || and § the expected numbers of signal and background events to be 
collected in a year. Using these numbers, statistical errors of the total cross section times 
branching fractions and the ratios of the branching fractions are calculated. The results are 
listed in Table [5]. As shown in this table, the results of the 2 jet analysis are comparable to 
those of the 4 jet analysis. The combined results show significant improvement compared 
to the results shown in Ref. It is found that the beam polarization gives rather modest 
improvement to the measurement of Br(H — > cc, gg) / Br(H — > bb). Similar studies are 
performed also for different Higgs masses. The statistical errors of Br(H — > cc, gg)/Br(H — ► 
bb) are shown in Fig. ^| as a function of the Higgs mass. At the low Higgs mass region, 
irreducible background from e + e~ — > Z°Z° increases, while the production cross section 
decreases in the high Higgs mass region. Thus the best result will be obtained if the Higgs 
mass lies around 110-120 GeV. 

Fig. 0A shows the relation between Br(H — > cc, gg)/Br(H — > bb) and m A o together with 
95 % confidence level lower limit for the SM Higgs boson with a mass of 120 GeV. If the 
ratio were measured to be 0.162 ± 0.028 with an integrated luminosity of 100 fb _1 (the SM 
prediction value with a relative error of 17 %), the MSSM which gives the ratio less than 
0.106 will be excluded at 95 % confidence level. In other words, a lower limit on m^o will be 
set to be 360 GeV at 95 % confidence level. To the contrary, if the ratio were measured to be 
significantly lower than the SM value, the SM will be excluded and both lower and upper 
limits on m A o will be obtained. Fig. shows such a case where the ratio were measured to 
be 0.100 with a similar relative accuracy. 
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Table 5: Expected measurement errors of various values after a year (50 fb x ) for the SM 
Higgs boson with a mass of 120 GeV (%). 
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Figure 7: Relation between Br(H — > cc, gg)/Br(H — > bb) and m A o for the Higgs boson with 
a mass of 120 GeV. (A) The solid (dotted) lines show 95 % confidence level lower limit after 
1 year (2 years). (B) The shaded area is the allowed region after a year (1 standard devia- 
tion). The region in the dotted lines shows the allowed region after two years. Kinematical 
thresholds for A direct production are also drawn. 
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Figure 8: The relation between Br(H — > cc, gg)/Br(H — > bb) and m&. The dotted 
lines show the uncertainty after two years. 

The fact that Br(H — > t + t~) can be measured with a statistical error of 13.4 % corre- 
sponds to that m b is determined with an error of 6.7 %. The systematic error of Br(H — > 
cc, gg)/Br(H — > 66) from the uncertainty of mt is estimated to be ~ 10 % after two years. 
FigH shows the relation between Br[H — > cc, gg)/Br(H — > 66) and m&. 

7 Summary 

We performed a detailed simulation study for the measurement of the ratio Br(H — > cc, gg) / Br(H 
bb) at -y/s = 300 GeV assuming the JLC standard detector. Here we introduce new ap- 
proaches to use all decay modes of Z° and to consider the beam polarization. We find that 
the results obtained from the 2 jet analysis are comparable to those of the 4 jet analysis, al- 
though the original event ratio is 2 jet : 4 jet ~ 20 : 70. Thus we can reduce the statistical 
error by combining the results of all decay modes. On the other hand, the effect of the beam 
polarization is found to be small for this analysis. 

The statistical error of the ratio Br(H — > cc, gg)/Br(H — ► 66) after two years is estimated 
to be 17 % for the SM Higgs with a mass of 120 GeV. With the good accuracy of the mea- 
surement, we have a chance to rule out the SM, or to constrain the mass of A in the MSSM, 
according to the measured value of the ratio. 

acknowledgments 

We would like to thank Prof. Yasuhiro Okada, Drs. Junichi Kamoshita, Minoru Tanaka and 
other members of the JLC Higgs Working Group for many useful discussions. We would 
like to also thank Drs. Junichi Kanzaki and Keisuke Fujii for great help for development of 
the detector simulator. 



13 



References 

[1] JLC Group, KEK Report 92-16, December, 1992. 

[2] C. Adolphsen, to be published in Proceedings of the Workshop on Physics and Experiments 
with Linear e + e~ Colliders, Morioka-Appi, Japan, September 8-12, 1995. 

[3] M. Leenen, to be published in Proceedings of the Workshop on Physics and Experiments with 
Linear e + e~ Colliders, Morioka-Appi, Japan, September 8-12, 1995. 

[4] J.P. Delahaye, to be published in Proceedings of the Workshop on Physics and Experiments 
with Linear e + e~ Colliders, Morioka-Appi, Japan, September 8-12, 1995. 

[5] V. Valakin, to be published in Proceedings of the Workshop on Physics and Experiments with 
Linear e + e~ Colliders, Morioka-Appi, Japan, September 8-12, 1995. 

[6] Y. Okada, M. Yamaguchi and T. Yanagida, Prog. Theor. Phys. 85 (1991) 1, 
Phys. Lett. B262 (1991) 54. 

[7] J. Ellis, G. Ridolfi and E Zwirner, Phys. Lett. B257 (1991) 83. 

[8] H.E. Haber and R. Hempfling, Phys. Rev. Lett. 66 (1991) 1815. 

[9] J. Kamoshita, Y. Okada and M. Tanaka, Phys. Lett. B328 (1994) 67. 
[10] J. Kamoshita, Y. Okada and M. Tanaka, KEK-TH-458, KEK Preprint 95-173. 
[11] M.D. Hildreth et al. , Phys. Rev. D49 3441 (1994). 
[12] Torbjorn Sjostrand, CERN-TH.7112/93 W5053/W5044. 

[13] T. Ishikawa, T. Kaneko, K. Kato, S. Kawabata, Y. Shimizu and H. Tanaka, "GRACE 
manual" KEK Report 92-19, February 1993. 

[14] JADE collaboration, W. Bartel et al. , Z. Phys. C33 (1986) 23. 



14 



